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Abstract 



We consider the production of the J/ip mesons in heavy ion colhsions at 
RHIC energies in the statistical coalescence model with an exact (canonical 
ensemble) charm conservation. The cc quark pairs are assumed to be created 
in the primary hard parton collisions, but the formation of the open and 
hidden charm particles takes place at the hadronization stage and follows the 
prescription of statistical mechanics. The dependence of the J /ip production 
on both the number of nucleon participants and the collision energy is studied. 
The model predicts the J /tp suppression for low energies, whereas at the 
highest RHIC energy the model reveals the J/ij: enhancement. 
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Typeset using REVTgX 



The experimental program for studies of the charmonium production in nucleus-nucleus 
(A+A) collisions at CERN SPS over the last 15 year was mainly motivated by a suggestion 
of Matsui and Satz to use the J/V' as a probe of the state of matter in the early stage 
of the collision. The original picture [|l| (see also Ref. for a modern review) assumes 
that charmonia are produced in primary collisions of the nucleons from the colliding nuclei. 
The number of created charmonium states is then reduced because of inelastic interactions 
with the nucleons of the colliding nuclei. An additional suppression may occur due to J/tp 
interaction with the secondary hadrons ('co-movers') 0. The probability to destroy the 
charmonium state increases obviously with the number of nucleon participants Np. Similar 
behavior is expected when the collision energy ^/s increases because the number of produced 
hadrons ('co-movers') becomes larger. This is known as the normal J/ip suppression. Fur- 
thermore, the J/ip (and other charmonia) are assumed to be formed mainly from the cc pairs 
with invariant mass below the D-meson threshold The fraction of these subthreshold 
pairs in the total number Nf^^ of cc pairs (which is roughly proportional to the number of 
produced open charm hadrons (D, D* , Ac etc.) also decreases with ^/s. Therefore, the ratio 

- ^ (1) 

is expected to decrease with increasing Np and/or ^/s. At large values of -y/i and Np the 
formation of the quark-gluon plasma (QGP) is expected which is supposed to be signaled 
by the anomalous suppression of the J/ip, i.e. a sudden and strong decrease of the ratio 
(0) is considered as a signal of the QGP formation. Hence, a decrease of the ratio (|l]) is an 
unambiguous consequence of the standard picture [0,0. 

A very different approach of the statistical J/ ip production, proposed in Ref. , assumes 
that J/ip mesons are created at the hadronization stage similar to other (lighter) hadrons. 

A picture of the J / ip creation via c and c coalescence (recombination) was subsequently 
developed within different model formulations [|6|-p!l|]. Similar to the suggestion of Ref. p, 
charmonium states are assumed to be created at the hadronization stage of the reaction, but 
they are formed due to the coalescence of c and c, which were produced by primary hard 
parton collisions at the initial stage. 

In this paper the Np and y/s dependences of the ratio (|l|) will be studied for Au+Au col- 
lisions at RHIC energies. We use the canonical ensemble (c.e.) formulation of the statistical 
coalescence model (SCM) []8|,p!I|. The number N^^^ of the produced cc pairs, which is the 
input for the SCM calculations, will be estimated within the perturbative QCD (pQCD). 
The considered pQCD-|-SCM approach reveals both the J/ip suppression (at N^c < 1) and 
the J/ip enhancement (at Ncc > 1) effects. 

In the framework of the ideal hadron gas (HG) model in the grand canonical ensemble 
(g.c.e.) formulation the hadron multiplicities are given by 



Irn] + k'^ — II j 
exp I ' ' ^ I ± 1 



(2) 
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where V and T correspond to the volume[] and temperature of the HG, rrij and dj denote 
particle masses and degeneracy factors. Eq.(^ describes the quantum HG: Bose and Fermi 
distributions for mesons and (anti)baryons, respectively. Quantum effects, however, are 
found to be noticeable for pions only, so that Eq.(|) for all other hadrons can be simplified 
to the Boltzmann result: 



where K2 is the modified Bessel function. 

In the case of complete chemical equilibrium the chemical potential /Xj in Eq.(^ is defined 

as 

/ij = bjfiB + qjfiQ + Sj^s + Cjfic , (4) 

where bj,qj, Sj,Cj denote the baryonic number, electric charge, strangeness and charm of 
hadron j. The baryonic chemical potential regulates the non-zero (positive) baryonic density 
of the HG system created in A+A collision. The chemical potentials fis aiid nc should 
be found as the functions of T and hb from the requirements of zero value for the total 
strangeness and charm in the system, and the chemical potential fiQ from the requirement 
of the fixed ratio of the electric charge to the baryonic number (this ratio is defined by the 
numbers of protons and neutrons in the colliding nuclei). 

The applications of the HG model to fitting the hadron abundances in particle and 
nuclear collisions revealed, however, a deviation of strange hadron multiplicities from the 
complete chemical equilibrium . It was suggested that strange quarks and antiquarks are 



distributed inside hadrons according to the laws of HG equilibrium, but the total number of 
strange quarks and antiquarks inside the hadrons is smaller than that in the equilibrium HG 
and remains approximately constant during the lifetime of the HG phase. Therefore, not only 
the "strange charge" Ns — N-s = but also the "total strangeness" Ng + Ng is then considered 
as a conserved quantity. In the language of thermodynamics, it means an introduction of an 
additional chemical potential, fi\s\, which regulates this new "conserved" number A^^ + A^. 
Then the additional term, (n^ + n^)fi\s\, should appear in the expression (^ for /ij , where 
ni and rz|- are the numbers of strange quarks and antiquarks inside hadron j. Introducing a 
notation, 7^ = exp{fi\s\/T) [|1^], one can implement this additional conservation according 
to the following simple rule: the hadron multiplicities A^ (|^) are multiplied by a factor 

^i^s+n-) ^ e.g. factor 7^ appears for K, K, A, A, S, S, factor 7^ for S, S and factor 7^ for fl, Q. 
For mesons with hidden strangeness, like rj, rj', uj, (f), having the wave function of the form 

Cu I uu) + Cd I dd) + Cs I ss) (5) 



To avoid complications we neglect the excluded volume corrections. The thermodynamical 
consistent way to treat the excluded volume effects was suggested in Ref. (see also [^] for 
further details) . If the excluded volume parameter is the same for all hadrons its effect is reduced 
only to the rescaling of the volume V: all particle number ratios remain the same as in the ideal 
hadron gas. 
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the factor 7^'*""'^ is used. 

From fitting the data on the hadron yield in particle and nuclear collisions it was found 
that 7s < 1 for all known cases. Parameter 7^ is called therefore the strangeness suppression 
factor. 

Recently an analogous procedure was suggested for charm hadrons |^. A new parameter 
7c has been introduced to treat simultaneously both the open and hidden charm particles 
within statistical mechanics HG formulation. The multiplicities Nj of single open charm 
and anticharm hadrons should be multiplied by the factor 7c and charmonium states by the 
factor 7^. In contrast to the suppression of strangeness in the HG (7^ < 1) one observes 
the enhancement of charm hadron yields in comparison to their equilibrium HG values. It 
means that 7c > 1 and this parameter is called the charm enhancement factor 

To take into account the requirement of zero "charm charge" of the HG in the exact 
form the c.e. formulation was suggested in Ref. In the c.e. formulation the charmonium 
multiplicities are still given by Eq.(|^) as charmonium states have zero charm charge. The 
multiplicities @ of the open charm hadrons will, however, be multiplied by the additional 



canonical suppression' factor (see e.g. |T^|). This suppression factor is the same for all 



individual single charm states. It leads to the total open charm multiplicity Nq in the c.e.: 

-o--o^,, (6) 

where Nq is the total g.c.e. multiplicity of all charm and anticharm mesons and (anti)baryons 
calculated with Eq.(^ and Io,h are the modified Bessel functions. For Nq « 1 one has 
Ii{No) / Io{No) — No/2 and, therefore, the c.e. total open charm multiphcity is strongly 
suppressed in comparison to the g.c.e. result. For Nq >> 1 one finds Ii{No) / Io{No) 1 
and therefore Nq Nq, i.e. the c.e. and the g.c.e results coincide. In high energy A+A 
collisions the total number of strange and antistrange hadron is much larger than unity. 
Hence the strangeness conservation can be considered in g.c.e. approach. The same is valid 
for baryonic number and electric charge. This is, however, not the case for the charm. At 
the SPS energies the c.e. suppression effects are always important: even in the most central 
Pb+Pb collisions the total number of open charm hadrons is expected to be smaller than one. 
It will be shown that the c.e. treatment of charm conservation remains crucially important 
also at the RHIC energies for the studies of the Np dependence of the open charm and 
charmonium production. Therefore, in what follows, the baryonic number, electric charge 
and strangeness of the HG system are treated according to the g.c.e. but charm is considered 
in the c.e. formulation where the exact "charm charge" conservation is imposed. 

Hence we formulate our model as follows. The charm quark-antiquark pairs are assumed 
to be created at the early stage of A+A collisions. In the subsequent evolution of the system, 
the number of this pairs remains approximately constant and is not necessary equal to its 
equilibrium value. The deviation from the chemical equilibrium should be taken into account 
by the charm enhancement factor 7c. The distribution of created cc pairs among open and 
hidden charm hadrons is regulated by statistical model according to Eq.(|^) with account for 
the canonical suppression (P). So the statistical coalescence model in the c.e. is formulated 
as: 

N"^' = - ^ No ^^^^^^^^ + 72 Nm (7) 
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where Nh is the total multiphcities of hadrons with hidden charm. (Note that the second 
term in the right-hand side of Eq.(^ gives only a tiny correction to the first term, i.e. most 
of the directly created cc pairs are transformed into the open charm hadrons.) To find Nq 
and Nh we use Eq.(|^) for the individual hadron thermal multiplicities in the g.c.e. and take 
the summation over all known particles and resonances [jl9[ with open and hidden charm, 
respectively^]. 

Provided that Nq, Nh and N^i'' are known, 7c can be found from Eq.(|^). The J/iIj 
multiplicity is then given by 

(^/^) = ll N'j% . (8) 

where iVj/*^ is given by 



NfU = Nj,^ + Bi{i,')N^, + Br(xi)iV^, + Br(x2)iV^, , (9) 



iVj/^, N^i, N^^, N^2 are calculated according to Eq.(||) and Br('?/'') ~ 0.54, Br(xi) — 0.27, 
Br(x2) — 0.14 are the decay branching ratios of the excited charmonium states into J/if'- 
Hence, to calculate the J /if) multiplicity in SCM we need the following information: 

1) the chemical freeze-out (or hadronization) parameters V,T,fiB for A+A collisions at high 
energies; 

2) the number N^^"^ of cc-pairs created in hard parton collisions at the early stage of A+A 
reaction. 

For the RHIC energies the chemical freeze-out temperature T is expected to be close to 
that for the SPS energies: T = 175 ± 10 MeV. To fix the unknown volume parameter V 
and baryonic chemical potential fiB use the parametrization of the total pion multiplicity 



where C = 1.46 GeV-^/^ ^nd rriTv is nucleon mass. For the RHIC energies the nucleon mass 
in Eq.(^) can be neglected so that 

(tt) CN, (v^)i/l (11) 

Eq.(PDD is an agreement with both the SPS data and the preliminary RHIC data in Au+Au 
collisions at i/i = 56 GeV and y/s = 130 GeV. The pion multiplicity ( p^Of ) should be equated 
to the total HG pion multiplicity A^^°* which includes the pions coming from the resonance 
decays (similar to Eq.(H)). The HG parameters V and are found then as the solution of 
the following coupled equations: 

(tt) =Ni°\V,T,f^B) ^ Vn':\T,f^B), (12) 
Np = VnB{T,fiB) , (13) 



^Note that possible (very small) contributions of particles with double open charm are neglected 
in Eq.(0). 
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where is the HG baryonic density. In these calculations we fix the temperature pa- 
rameter T. The baryonic chemical potential for Au+Au collisions at the RHIC energies is 
small {fiB < T) and decreases with collision energy. Therefore, most of the thermal HG 
multiplicities become close to their limiting values at /i^ — >^ 0. Consequently most of hadron 
ratios Nj°^/Nj^°^ become independent of the collision energy. The volume of the system is 
approximately proportional to the number of pions: 



Note that T = 170 180 MeV leads to the HG value of the thermal ratio: 



3/2 



exp 



T 



0.04^0.05 



(14) 



(15) 



in agreement with data [0 in central {Np > 100) Pb+Pb collisions at the CERN SPS. This 
fact was first noticed in Ref. []T^. At small Np as well as in p+p and p+A collisions the 



measured value of the {4'')/ {J/'ip) ratio is several times larger than its statistical mechanics 
estimate (^). Our analysis of the SCM will be therefore restricted to A+B collisions with 
Np > 100. We do not intend to describe the open and hidden charm production in p+p, 
p+A and very peripheral A+B collisions within the SCM. 

The number of directly produced cc pairs, N^^"^, in the left-hand side of Eq.(0) will be 
estimated in the Glauber approach. For A + B collision at the impact parameter b, this 
number is given by the formula: 



Ar^^-(6) = AB TAsip) a^NN -^cc + X) 



(16) 



where a{NN cc + X) is the cross section of cc pair production in + collisions and 
Tab{P) is the nuclear overlap function (see Appendix for details). 

The cross section of cc pair production in N+N collisions can be calculated in the pQCD. 
Such calculations (in the leading order of the pQCD) were first done in Ref. |2^. We 
use the next-to-leading order result presented in Ref. [^]. This result was obtained with 
GRV HO structure functions, the c-quark mass and renormalization scale were fixed 
at TTif. = /X = 1.3 GeV to fit the available experimental data. We parametrize the ^/s- 
dependence of the cross section for = 20 -i- 200 GeV as: 



a{pp 



cc) 



(To 1 




(17) 



with (Jo ^ 3.392 fjh, Mq 2.984 GeV, a ^ 8.185 and p ^ 1.132. Note that free param- 
eters of the pQCD calculations in Ref. |2^ were fitted to the existing data, therefore, our 
parametrization {\L7\) is also in agreement with data on the total charm production in p+p 
collisions. 

The average number of participants ('wounded nucleons') in A+B collisions at impact 
parameter b is given by 



Np{b) = A 

+ B 



+ 00 /• + 00 , 

dx / dy TAiyx^ + y^] 



dx 

OO J — OO 



dyTBUx^ + {y-bY 



l-[l-aT^TB{^x^ + {y-bf) 
1 - ( 1 - a^NTAiJ^y^ 



;i8) 
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where Ta{s) (Tb{s)) is the nuclear thickness function for the nucleus A (B) and a]^^ is the 
total inelastic cross section of + interaction. To parametrize the i/s-dependence of a 



inel 
NN 

we made the assumption that in the energy range ^/s = 10 -h 200 GeV, it is proportional to 
the total NN cross section aNN and use the standard fit for (Xnn 

cf^t - O.TorjvTv ^ 0.7 (Xpps^ + Yvp.s-'^' - Y2vvs"'''). (19) 
where e = 0.093, r/i = 0.358, r^a = 0.560, Xpp = 18.751, Y^y^ = 63.58 and Y2pp = 35.46. 



Eqs. (|lg) and (|18D provide parametric dependence of the number of produced cc pairs 
on the number of participating nucleons, N^-^'^ = N^^^{Np), which is shown in Fig |I| for 
Au+Au collisions. It is seen that the dependence is represented by a straight line in the 
double-logarithmic scale, so that 



N^' ~ (Np)' (20) 

for Np > 50. We find that k = 1.31 ^ 1.350. Using Eq.(|I^ one finds then the following 
behavior of N^^'^ at high energies: 

- (iVp)'(v^)" • (21) 

Now we are able to calculate the ratio R (|1|) in the SCM and study its dependence on 
Np and y/s. The dependence of R on the number of participants is shown in Fig.|. It is 
seen that the ratio has qualitatively different behavior at different energies. At the lowest 
RHIC energy ^/s = 56 GeV, the SCM predicts decreasing of the ratio with the number of 
participants {J/i' suppression). In contrast, at the highest RHIC energy y/s = 200 GeV 
the ratio increases with the number of participant {J/ip enhancement) for Np > 100. Both 
suppression (at Np < 150) and enhancement (at Np > 200) are seen at the intermediate 
RHIC energy = 130 GeV. 

Similarly, there are qualitatively different dependencies of R on the collision energy for 
small {Np = 100) and for large {Np = 350) number of the participants. This can be seen in 
Fig.^. Non-monotonic dependence of the ratio R on ^/s is expected at Np = 100 . At A'p = 
350, the ratio R increases monotonically with ^/s at all RHIC energies ^/s = 56 -h 200 GeV. 
The minimum of R in this case corresponds to the energy region between the SPS and RHIC: 

^ 30 GeV. 

To understand the behavior of R it is instructive to study the limiting cases: N^-^"^ « 1 
and N^^"^ >> 1. Neglecting the hidden-charm term in Eq.(|^) one finds for N^^^ << 1 : 



hence, 



Nf ^ 1 N'o , (22) 



N^^ N^/A V 



^It is interesting to note that for the most central ^ + yl collisions {Np ~ 2 A) , N^^^ has approxi- 
mately the same dependence on the atomic weig ht of the colliding nuclei: iVj*^ ~ A'^l'^ ~ {Npf/^ ■ 
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Eq.(p3D shows XjV universal suppression of the ratio R. This ratio decreases as A^^ ^ and 



{^/s) with the increasing number of participants and colhsion energy, and the shape of 
this J/ip suppression is essentially independent of the functional dependence of N^-^^ on Np 
and ^/s. 

If A^J*" » 1 one finds from Eq.(0): 

^ lie No, (24) 
so that 7c ~ 2N^I'/No ~ N^^' /V and, hence, 

R^mi ^ MSi ^ M; ^ N"-^ i^)'"" (25) 



According to Eq. (p5|) the ratio R increases with both Np and y/s. The J/ip enhancement 
takes place due to the fact that the number of primary nucleon-nucleon collisions grows faster 
than the number of participants (A^^*'' ~ (Np)'', k > 1) and because the pion multiplicity (and 
therefore the volume of the system) is less sensitive to the collision energy ((tt) ~ (v^)^^^) 
than the number of cc pairs (A'^*'' ~ (-\/i)^, P > 1/2). 

It is seen from Fig.|I| that A"^"' << 1 at the lowest RHIC energy for small numbers 
of participants, hence the SCM predicts the J/ip suppression. In contrast, for the highest 
RHIC energy and large Ap the opposite limit A"^*^ >> 1 is reached. This leads to the J/^jj 
enhancement. 



In conclusion, the production of the J/ip mesons is studied in Au+Au collisions at 
the RHIC energies in the statistical coalescence model with the exact charm conservation. 
The cc quark pairs are assumed to be created in the primary hard parton collisions and 
their number is estimated within the pQCD. At the hadronization stage the cc quarks are 
distributed among the open charm and charmonium particles according to the hadron gas 
statistical mechanics in the canonical ensemble formulation. 

Decreasing of the (J/ip) to N^^"^ ratio with increasing the number of nucleon participants 
Np is found at the lowest y/s = 56 GeV RHIC energy. At fixed small number of partici- 
pants {Np ^ 100) the ratio decreases with a/s between the lowest (i/i = 56 GeV) and the 
intermediate {^/s = 130 GeV) RHIC energies. This is in a qualitative agreement with the 
standard picture [jl],^ of the J/ip suppression. In contrast, a rise of the {J ftp) to N;!!'' ratio 
with the collision energy is predicted for central Au + Au collisions. Moreover, at the highest 
RHIC energy, the ratio is expected to grow with the number of participants, which is in a 
drastic contradiction with the standard picture. The reason for this that in the standard 
picture the hidden charm mesons are supposed to be created exclusively in the primary 
(hard) nucleon-nucleon collisions. It is assumed that all other interaction can only destroy 
them. Especially strong suppression of the charmonia is expected in the quark-gluon plasma 
('anomalous J/ip suppression'). In distinction to this standard approach, the statistical co- 
alescence model considers a possibility for the charmonium states to be formed from c and c 
at the stage of the quark-gluon plasma hadronization. This possibility definitely cannot be 
ignored, when the number of produced cc pairs per A + A collision becomes large: Ncc » I 
(this happens for the central Au + Au collisions at the highest RHIC energy). In this case 
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the c and c pairs initially produced in different hard collision processes can recombine into a 
hidden charm meson. Therefore, an increase of the (J/ip) to N^}^"^ ratio should be expected. 

The hot quark gluon plasma is most probably formed at high RHIC energies and this 
destroy all primarily produced charmonium states [^]. However, the hadronization of the 



quark gluon plasma within the SCM reveals itself in the J/ip enhancement rather than 
suppression. Another interesting phenomena may also take place: when the number of cc 
pairs becomes large, two c quarks (or two c) can combine with a light (anti) quark and form 
a double charmed (anti)baryon. These baryons are predicted by the quark model but have 
not been observed yet. We expect that double (and probably triple) charmed baryons may 
be discovered in the Au + Au collisions at RHIC ||2l 
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APPENDIX: NUCLEAR GEOMETRY 

The spherically symmetrical distribution of the nucleons in the Au — 197 nucleus can be 
parametrized by the two-parameter Fermi function p8[ (this parametrization is also known 
as Woods-Saxon distribution): 



p{r) = Po 



, r — c^'' ^ 

1 + exp 



(Al) 



a 

with c ~ 6.38 fm, a ~ 0.535 fm and po is given by the normalization condition: 

Air / drr'^pir) = 1. (A2) 



The nuclear thickness distribution Ta(6) is defined by the formula 



TAib) = / dzp(V¥+^) , (A3) 
and the nuclear overlap function is defined as 

TABib) = J^^ dx dyTA ijx^ + y^) Tb [^Jx^ + {y~hY^ . (A4) 



From Eq.(|A^), one can deduce that the above functions satisfy the following normalization 
conditions: 

/•oo roc 

271 dbbTA{b) = l, 2-n dbbTAB{b) = l. (A5) 
Jo Jo 
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FIG. 2. (J/V')/A^J'' versus Np for = 56,130,200 GeV. The vertical line shows the lower 
bound of the applicability domain of the SCM. 
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FIG. 3. {J/'^)/N^f versus ^/s for Np = 100 and 350. 
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